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Stimulated by the anomalous superconducting properties of UPts, we investigate the 
pairing symmetry and the transition temperature in the two-dimensional(2D) and three- 
dimensional(3D) hexagonal Hubbard model. We solve the Eliashberg equation using the 
third order perturbation theory with respect to the on-site repulsion U. As results of the 
2D calculation, we obtain distinct two types of stable spin-triplet pairing states. One is the 
/-wave(Bi) pairing around n — 1.2 and in a small U region, which is caused by the ferro- 
magnetic fluctuation. Then, the other is the Px{ot: Pj^)-wave(Ei) pairing in large U region 
far from the half- filling [n = 1) which is caused by the vertex corrections only. However, 
we find that the former f-wave pairing is destroyed by introduced 3D dispersion. This is 
because the 3D dispersion breaks the favorable structures for the /-wave pairing such as 
the van Hove singularities and the small pocket structures. Thus, we conclude that the fer- 
romagnetic fluctuation mediated spin-triplet state can not explain the superconductivity of 
UPta. We also study the case of the pairing symmetry with a polar gap. This p2-wave(Ai) 
is stabilized by the large hopping integral along c-axis tz- It is nearly degenerate with the 
suppressed (or pj,)-wave(Ei) in the best fitting parameter region to UPta (1.3 <tz< 1.5). 
These two p-wave pairing states exist in the region far from the half-filling, in which the 
vertex correction terms play crucial roles like the case in Sr2Ru04. 

KEYWORDS: UPta, third-order perturbation theory, spin-triplet superconductivity, heavy 
fermion 



1. Introduction 

Much attention have been focused on hexagonal heavy fermion superconductor UPta, 
since it was discovered in 1984.^ The heavy fermion state is characterized by its large value 
of the electronic specific heat coefficient 7 = 420mJ/K^moleU. The coefficient of T'^-law of 
resistivity is also enhanced by 10^ times as large as that of conventional metal. Such behavior 
arises from the strong electron correlation. In UPta, the itinerant electrons compose such 
a Fermi liquid state and then undergo superconducting transition. In zero field, the two 
superconducting transitions occur at ~0.48 K and ~0.52 K'^.^ The field-temperature phase 
diagram was constructed experimentally'' and it consists of at least three phases (denoted as 
A, B and C). Most of theoretical models for this splitting of Tc are based on the coupling 
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between the superconducting order parameter and the smaU antiferromagnetic moment (us = 
0.02^5 /U) which appears below Tn 5K^.® 

Furthermore, a power law behavior of the ultrasound attenuation,^ NMR relaxation rate,^ 
specific hcat^'^ and penetration depth^ below Tc suggested an unusual order parameter with 
energy gap which vanishes at some points on the Fermi surface. In addition, the NMR^^*^ and 
muon spin rotation-relaxation (^uSR)^^ experiments detected no change of the Knight shift 
below Tc. These suggest a possibility of the odd-parity of the Cooper pairing. 

In order to explain these unconventional properties in UPts, various phenomenological 
theories have been proposed. Theories based on two-dimensional representation-*^^ -"^'^ assign a 
gap with both point nodes along c-axis and a line of nodes at = (called a hybrid gap) 
to the phase B. On the other hand, the theory based on one-dimensional representation^^ 
requires either a hybrid gap or a gap with only a line of nodes at fcz = (a polar gap). The 
experiments of thermal conductivity^^ and penetration depth^ suggested that a hybrid gap 
is favored, but the gap structure of UPts has not been completely determined. 

In contrast with the phenomenological theories, few microscopical theories determin- 
ing the pairing symmetry and the gap structure have been reported. The band calcula- 
tions"'^^^^,"'^^ based on the assumption that 5f-electrons are itinerant, reported that five three- 
dimensional bands form the Fermi surface. The results of the de Haas-van Alphen (dHvA) 
experiments^^ ^^^'^ are consistent with these band calculations, and the observed effective mass 
of each band is almost 20 times as large as the corresponding calculated mass. Among these 
five bands, we take up band 37-clectron in rcf.l9 (a rugby ball like band structure with small 
electron pockets at K-points in the hexagonal Brillouin zone) because it has the heaviest effec- 
tive mass and occupies the largest area of the Fermi surface. Then, in the following sections, 
we discuss the pairing symmetry in the microscopic way based on band 37-electron. 



2. Model 



We consider the following single band Hubbard model on hexagonal lattice, 



H 



Ho + ilint 



(1) 




(2) 



ki 

From the tight binding approximation, energy dispersion ^(/c) is given by 




(3) 
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(4) 



Here, ti, t2 and ^3 are the nearest-, the next-nearest-, and the third- nearest-neighbor hopping 
integrals in ab- plane respectively. The tz is the hopping integral along c-axis, and /x is the 
chemical potential. Now we fix ti = 1, and calculate physical quantities changing the other 
variables t2, ^3, tz-, U and electron filling n. 

3. Formulation 

One of the authors (K.Y.) has described the Fermi liquid theory for the heavy electron sys- 
tem^^.^^ Based on this idea, we start from the renormalized quasiparticle state, and calculate 
the momentum dependence of the effective interaction between them. 

Recently, some authors have calculated the transition temperture Tc and the pairing sym- 
metry for many superconductors including high-Tj. cuprates and Sr2Ru04 on the basis of the 
third order perturbation theory (TOPT) with respect to U?'^ All of them are in good agree- 
ment with experimental facts. Although the convergence of the vertex correction terms is an 
important problem, TOPT in the moderately correlated region has been justified by higher 
order calculations.^^ Thus, in the same framework, we apply TOPT to our model. 

In TOPT, the effective interaction for the spin-singlet state is given by 



v^{k, k') = v|pA(fc, k') + y4tex(fc, k') 

U + U\oik-k') + 2U^xlik-k') 
2C/3Re G{k + q)G{k' + q)[xo{q) - Mq)] 



^rpa(^5 k' 



Vertex (^5 k') 



For the spin-triplet state, 

V^{k,k') 



^RPa(^5 k') + K^rtex(^) 

V^Mk,k') = -U\o{k-k') 

2U^Re[j2 G{k + q)G{k' + q)[xo{q) + Mq)] 



Vertex 



{k,k') 



where 



Xo{q) 

G{k) 



Y,G{k + q)G{k) 
k 

Y,G{q-k)G{k) 



(5) 
(6) 
(7) 

(8) 
(9) 
(10) 

(11) 
(12) 

(13) 



iun - i{k) 

Here, is a shorthand notation as A; = (fc, iun) and ojn = 7rT(2n -|- 1) is a fermion Matsubara 
frequency, ^(k) is the band of quasiparticles with large effective mass. In Fig.l, the diagrams 
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Fig. 1. Feynman diagrams of the effective interaction up to third order. Solid and dashed hnes cor- 
respond to the Green's function and the interaction, respectively. 



of the effective interaction are shown. The diagrams enclosed by a dashed line are the ver- 
tex correction terms. We calculate Tc and the pairing symmetry by solving the Eliashberg 
equation, 



X^^A{k) = --^Vik,k')\G{k')\''A{k'). 



(14) 



At T = Tc, the maximum eigenvalue Amax becomes unity, and the pairing symmetry is deter- 
mined by the momentum dependence of the gap function A(A;). Thus, by estimating Amax) we 
can determine which type of pairing symmetry is stable. Here, we have ignored the damping 
effect due to the normal selfenergy. This contribution is important for estimating Tc quanti- 
tatively, but qualitative properties of superconductivity such as the pairing symmetry is not 
usually influenced. 

4. Results of 2D Calculations 

In this section, we show the results of two-dimensional calculations, considering the kz = 
plane (including F, K, M points in the hexagonal Brillouin zone) of the band 37, that is, we 
fix tz = 0. Here, we notice that band 37 is highly three-dimensional, but the calculations on 
the basal triangular lattice gives us much information as starting points. We will show the 
results of 3D calculations by introducing tz in §5. Then, we also fix = and T = 0.01 in this 
section. The calculated Fermi surfaces are shown in Fig.2. We can see larger electron pockets 



4/15 



J. Phys. Soc. Jpn. 



Full Paper 



(a) (b) 
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n=0. 8 

n=1.0 
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Fig. 2. Calculated Fermi surfaces in the case of (a)t2 = 0.3, n ~ 0.6,0.8,1.0,1.2 (b)t2 = 0.4, n = 
0.6,0.8,1.0,1.2 {c)t2 = 0.5, n = 0.6,0.8,1.0 {d)t2 = 0.6, n = 0.6,0.8,1.0 . Electron pockets at 
K-points become large as either n or t2 increases. 

at K-points as either n or t2 increases. Within the kz = plane, the best fitting parameter 
region for UPts is t2 = 0.4 ~ 0.6, n = 0.7 ~ 1.0. The possible candidate of pairing symmetry 
in 2D triangular lattice belongs to the following irreducible representations of Cq^. 



symmetry One of basis function 



Bi (/i-wave) sin ^ [cos - cos 
B2 (/2-wave) sin ^ (cos ^ - cos ^) 



El (p- wave) 



sin cos ^ 
cos sin ^ + sin ky 



, sm^sin-2 
E2((i-wave) \ /^i. k 

' cos ^^^^ cos ^ - cos h. 



V 



First, we set [/ = 2. Figure 3 is a contour plot of Amax as a function of t2 and electron 
filling n (n = 1 at the half filling). The black thick line represents the crossover line of two 
pairing symmetries. When n is larger than half filling, an /i(Bi)-wave spin-triplet state is 
stable. This /i-wave state was predicted by Ikeda et al?^ in the case of NaxCo02 • yH20. 
In contrast, for the case less than half filling, d(E2)-wave pairing state covers a wide region. 
But the eigenvalues of this d-wave state is too low to be realized. Then, what is the origin of 
/-wave pairing above half filling? In figure 4, we show the bare spin susceptibility xo (9*, 0) in 
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Fig. 3. The contour plot of Amax as a function of t2 (band width) and n in the case of f = 2 and 
T = 0.01. An /i-wave spin-triplet state is stable in the large n region. The eigenvalue for rf-wave 
pairing below the half filling is too low to be realized. 

the case of t2 = 0.3, n = 1.2. The peak structures can be seen at F-points, and cannot be seen 
in the region where d-wave pairing state is stable. This peak corresponds to the increase of 
the density of states by van Hove singularity, which can be seen around K-point in Fig.2(a) 
(n=1.2). The peak structure at F-point implies the ferromagnetic instability. These conditions 
stabilize the /i-wave triplet pairing. Next, to examine influence of the vertex corrections, we 
calculate the eigenvalues by including only the RPA-like diagrams of the effective interaction 
up to the third order. The n-dependence of the calculated eigenvalues for /i-wave pairing in 
the RPA-like calculation, and that for /i- and d-wave in TOPT are shown in Fig.5. We can 
see the vertex correction terms work in favor of realizing the triplet pairing, as Nomura et al. 
predicted in Sr2Ru04.^^ Here, we have to note that the Fermi surface in this /i-wave region 
is rather larger than actual one. This does not continue to the best fitting region for UPta 
{t2 = 0.4 ~ 0.6, n = 0.7 ~ 1.0), because of the absence of the van Hove singularities near the 
Fermi surface in spite of the existence of the electron pockets. 

Next, we set U = 5 to investigate the triplet pairing in the best fitting parameter region. 
Figure 6 is a contour plot of Amax as a function of t2 and n. We can see p(Ei)-wave spin-triplet 
pairing state appears in low n region. In this region, there are no significant peak structures at 
F-point, because van Hove singularities and/or small pocket structures are lacking. Thus, this 
p-wave pairing is caused by the vertex correction only, in contrast with the previous /-wave 
pairing. Then we increase U ioU = 7, and show a contour plot of Amax as a function of t2 and 
n in Fig. 7. We can see the p-wave pairing becomes stable. Far from the half-filling, the p-wave 
pairing covers a wide region which includes the best fitting parameter region for UPt3. Here, 
we note, in the case of n > 0.9 for large U, we cannot perform reliable calculations because 
U maybe too large compared with 1/x™^, where 1/x™^ is the maximum value of the bare 
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Fig. 4. The bare spin susceptibility xo(9, 0) in the case of t2 = 0.3, n = 1.2 and T = 0.01. We can 
see peak structures around F-points. 
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Fig. 5. The n-dependence of eigenvalues in the case of t2 = 0.3, U = 2 and T = 0.01. The dashed 
(dotted) line with black (white) squares is the result for /i-wave (d-wave) pairing in TOPT. The 
line with black circles is the result for /i-wave pairing in the RPA-like calculation. 



susceptibility Xo(^) 0). 

In this section, we have obtained distinct two types of stable triplet pairing states. One is 
the /(Bi)-wave pairing in large n and small U region, which is caused by the ferromagnetic 
fluctuation. The vertex correction terms also increase Amax for this /-wave pairing. Another 
is the p(Ei)-wave pairing in large U region far from the half filling. This state is caused by 
the vertex corrections only, but the form of the Fermi surface is in better agreement with that 
from the band calculations. However, since the band 37 is highly three-dimensional, we have 
to introduce 3D dispersion and investigate the stability for the effect of three-dimensionality. 
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Fig. 6. The contour plot of Amax as a function of t2 (band width) and n in the case oi U = 5 and 
T = 0.01. Here, we consider the case of less than half-filling only. A p-wave pairing appears far 
from the half-filling. 




Fig. 7. The contour plot of Amax as a function of t2 (band width) and n in the case oi U = 7 and 
T = 0.01. A p-wave triplet pairing is stabilized by the vertex corrections for large U, and covers 
the wide region. 



5. Results of 3D Calculations 

5. 1 Stability for 3D dispersion 

In this section, we investivate the effect of 3D dispersion on these two pairing states, 
introducing tz- We also &x ts = and T = 0.01 in this section. First, we apply the 3D 
calculation to the high-filling /i(Bi)-wave pairing state at U = 2. In Fig.8, we show the in- 
dependence of eigenvalues for /i-wave state in two cases t2 = 0.3, n = 1.2 and t2 = 0.3, n = 
1.15. As tz increases, the eigenvalues show relatively sharp decrease (compared with the next 
p-wave cases) and then vanishes near tz = 0.3. This result is different from the case of d- 
wave pairing on the square lattice:^^ In that case Tc begins to be suppressed at tz = 0.4, 
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Fig. 8. The t^-dependence of eigenvalues for /i(Bi)-wave pairing in the case of t2 = 0.3, U = 2 and 
T = 0.01. The Hne with black squares (circles) is the result for n = 1.2 {n = 1.15) in TOPT. The 
suppression of Amax by 3D dispersion is very strong. 



then decreases by one order of magnitude at tz = 1.0. In our case, the rapid decrease of 
the eigenvalue for /i-wave indicates that the ferromagnetic fluctuation mediated spin-triplet 
superconductivity is easily destroyed by three-dimensionality, which breaks both the van Hove 
singularities and the small pocket structures. Actually, the band 37 is highly three-dimensional 
{tz > 1) and both structures exist only near the kz = plane. Thus, we conclude that this 
/i-wave pairing can not explain the superconductivity of UPts. 

Next, we apply the 3D dispersion to the low-filling p(Ei)-wave state aiU = 7. In Fig.9, we 
show the t^;-dependence of eigenvalues for p-wave state in the two cases t2 = 0.4, n = 0.6 and 
t2 = 0.4, n = 0.7. In contrast to the /i-wave pairing case, the decrease of the eigenvalue of the 
p-wave pairing is rather slow, especially far from the half-filling. This is because the p-wave 
pairing is not caused by the special structures of Fermi surface but the vertex corrections. 
When we introduce tz, the calculated band width become larger. Thus we increase the on-site 
repulsion U and try to obtain the reliable eigenvalues in the next subsection. 

5.2 Pairing Symmetry with a Polar Gap 

There are many pairing symmetries in hexagonal lattice {Dq^) in addition to the irreducible 
representations of Ce^, for example ^ig{dxz and dyz are degenerated), ^2^{dxy and d^2_y2 are 
degenerated). Big (5- wave), B2g(5-wave) for spin-singlet state, Ai^{pz), Eiu(pj; and py are 
degenerated), E2n{fxyz and j(x2-y2)z are degenerated), Biu(/i), B2u(/2) for spin-triplet state. 

In this section, we investigate Amax using 3D band structure from beginning. Calculated 
Fermi surface in the case of t2 = 0.3, = —0.5, n = 0.45 and tz = 1.4: is shown in Fig. 10. This 
well describes the band 37-electron. Such a low electron filling (n = 0.45) is due to the absence 
of electron in kz = ±7r plane. In Fig.ll, we show the ?7-dependence of positive eigenvalues at 
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Fig. 9. The ^2-dependence of eigenvalues for p(Ei)-wave pairing in the case of t2 = 0.4, U = 7 and 
T = 0.01. The soUd (dashed) hne with black squares (circles) is the result for n = 0.6 (n = 0.7) in 
TOPT. This low-filling p-wave pairing is stable for 3D dispersion compared with the /i-wave in 
Fig.8. 



T = 0.01. Because we consider the region far from the half-filling, the d-wave pairing states 
are suppressed by the large vertex corrections from large U. On the other hand, the p-wavc 
pairings are stabilized by them. The two p-wave pairing (Aiu and Eiu) are nearly degenerate, 
and the stable state between them is changed by the small variations of some other parameters 
such as t2, t^, n and t^. Thus it is difficult for us to determine which p-wave pairing is stable. 
Here, we note that the band width in this parameter equals 16.0, which is much larger than 
that in 2D calculations, thus such large value of U in Fig. 11 is usable. Next, in Fig. 12, we show 
the t^-dependence of eigenvalues for p2:(Eiu)-wave and p2(Aiu)-wave pairing states at U = 9 
and T = 0.01. The eigenvalue for pa;(Eiu)-wave pairing is suppressed by 3D dispersion as 
shown in §5.1. On the other hand, p2(Aiu)-wave pairing is stabilized by tz- These two p-wave 
states are nearly degenerate in the best fitting parameter region to UPts (1.3 < tz < 1.5). 
The pziAi u)-wave pairing agrees with experimental results which confirm the gap structure 
with a line of nodes at kz = plane. However, we cannot explain the E2u hybrid gap (fxyz 
and f(x^-y2)z)- We may have to use more realistic models such as a multi-orbital model to 
explain it. 

6. Discussions 

In this paper, we have investigated the single band Hubbard model for UPta on the 
basis of the perturbation theory. We have performed all the calculations in renormalized 
scheme, where we consider ^(k) to be the dispersion of quasiparticles with heavy effective 
mass. In that scheme, we should not include the normal selfenergy corrections giving rise to 
mass enhancement which has been already included by the renormalization. Therefore we 
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Fig. 10. The Fermi surface at t2 = 0.3, ts = —0.5, n = 0.45 and tz = 1.4 (a)in ky — plane, (b)in 
kz = plane. This well describes the band 37-electron in ref.l9. 




Fig. 11. The iJ-dependence of eigenvalues for p^-wave (Aiu) pairing and Px{py)-waye (Eiu) pairing 
in the case of t2 = 0.3, = -0.5, n = 0.45, tz = 1.4 and T = 0.01. The band width is 16.0. The 
solid (dashed) line with black squares (circles) is the result for px-wave (pz-wave) in TOPT, and 
the dotted line with black triangles is the result for dxyid^^ _y2)--wavc (E2g). In large U region, 
two p-wave pairings are stabilized by the vertex corrections, and these are nearly degenerate. 



have ignored them. As a result, we have also ignored the damping effect due to the normal 
selfenergy, but qualitative properties of superconductivity such as the pairing symmetry is 
not usually influenced. 

When we examine T(. quantitatively, we have to estimate the band width of the heavy 
quasiparticles. This is very difficult especially in the case of the heavy fermion systems. The 
effective mass and the band width contribute directly to the value of Tc, as Sasaki et al. 
account for the difference of Tc for high-Tp cuprates of La system and that of Y system 
using the d — p model.^^ Thus we have to adopt more realistic models such as the periodic 
Anderson model and calculate the higher order selfenergy corrections, to examine both the 



11/15 



J. Phys. Soc. Jpn. 



Full Paper 



1.0 



...... 



Pz(Aiu)-wave - 



A 



0.2 



0.8 



0.6 



0.4 





0*^ 



(10. 4) 



0.2 0.4 0.6 0.8 



1 

(14.4) 



12 14 1 .6 
(16.8) 



tz 

(band width) 



Fig. 12. The i^-dependence of eigenvalues for p^-wave (Aiu) pairing and Pa;-wave (Eiu) pairing in 
the case of t2 = 0.3, is = -0.5,n = 0.45, J7 = 9 and T = 0.01. The soUd (dashed) line with black 
squares (circles) is the result for p^^-wave (px-wave) in TOPT. The Pz-wave pairing stabilized by 
tz is nearly degenerate with the suppressed Px(or Py)-wave in the region 1.3 <tz< 1.5. 

heavy effective mass and Tc quantitatively. 

Yamagami also calculated the energy band structure for UPts^^ using a fuUy-relativistic 
spin-polarized version of the linearized augmented-plane-wave (LAPW) method. The calcu- 
lated energy band structure show us fact that, compered with the one we calculated in sub- 
section 5.2, the six pockets are large and the main rugby ball like Fermi surface locates on the 
fiat bottom of the dispersion. Thus, we have introduced further neighbor hopping integrals to 
reproduce his energy band structure and then calculated eigenvalues more precisely. However, 
the calculated density of state at the Fermi surface and the calculated Tc do not distinctly 
differ from those in subsection 5.2. Therfore, we do not miss the essence of band 37-electron 
using the simple Fermi surface calculated in subsection 5.2. 

Our simple calculations predict the following points. 

(1) The d-wave pairing is not possible because the antiferromagnetic fluctuation is suppressed. 

(2) Only the triplt (odd parity) pairing is probable. 

(3) Simple ferromagnetic fluctuation machanism can not explain the superconductivity of 
UPt3. This is because the van Hove singularities and the small pocket structures exist 
only near the kz = plane. 

(4) The p-wave pairing is the most probable. It is stabilized by the large vertex corrections 
from large U. 

(5) The calculated Tc for p-wave pairing is rather low. In the renormalized scheme, the limit of 
the effective on-site repulsion is not clear because of the ambiguity of the band width of the 
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quasiparticles. In some case, we need to use larger value to realize the superconductivity 
of UPtg. 

7. Conclusion 

In conclusion, on the basis of the microscopic theory, we have investigated the pairing sym- 
metry and the transition temperture in the two-dimensional(2D) and three-dimensional(3D) 
hexagonal Hubbard model considering the band 37-electron in UPta. We have solved the 
Eliashberg equation using the third order perturbation theory with respect to U. As the re- 
sults of the 2D calculation, we obtain distinct two types of stable spin-triplet pairing states. 
One is the /-wave(Bi) pairing in large filling n and small U region, which is caused by the 
ferromagnetic fluctuation. Then the other is the Px{ov P2^)-wave(Ei) pairing in large U region 
far from the half filling which is caused by the vertex corrections only. However, we find that 
the former /-wave pairing is easily destroyed by introduced 3D dispersion. This is because 
the 3D dispersion breaks the favorable structures for the /-wave pairing such as the van Hove 
singularities or the small pocket structures. Thus, the ferromagnetic fluctuation mediated 
spin-triplet state can not explain the superconductivity of UPta. We have also studied the 
case of pairing symmetry with a polar gap. This p2-wave(Ai) is stabilized by both the large 
hopping integral along c-axis tz and large U. It is nearly degenerate with the suppressed Px{oT 
Py)-wave(Ei) in the best fitting parameter region to UPts (1.3 < tz < 1.5). The Pz-wave 
pairing partially agree with the experimental results. Finally, we note that these two ^?-wave 
pairing states exist far from the half-filling, in which the vertex correction terms play crucial 
roles like the case in Sr2Ru04. 
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